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Abstract: A study of two-Higgs-doublet models (2HDM) in the decoupling limit reveals 
the existence of parameter configurations with a large triple-Higgs self-coupling as the 
only low-energy trace of the departure from a Standard Model (SM) Higgs sector. This 
observation encourages attempts to search for double Higgs production at the Large Hadron 
Collider (LHC) even in mass regions which have been shown to be very hard to probe in 
the context of SM-like Higgs self-couplings. In this document we consider the case of an 
Intermediate Mass Higgs (IMH) boson, with 120 GeV < Mh < 140 GeV, produced in 
pairs via vector-vector fusion, Higgs-strahlung and associate production with heavy-quarks 
and decaying into bb pairs. After a detailed signal-to-background analysis, we confirm that 
the observation of a Higgs-pair signal is very challenging in the framework of the SM, even 
at the LHC with upgraded luminosity. In contrast, we verify that the sensitivity is sufficient 
to detect departures from the SM which would escape detection in the measurements of 
the single-Higgs production channels. 

Keywords: Electro- Weak Symmetry Breaking, Higgs potential, Hadron Colliders, 
Beyond the Standard Model. 



*On leave of absence from: Dipartimento di Fisica Teorica, Torino and INFN Sezione di Torino, Italy 
^Present address: Centro Studi e Ricerche "E. Fermi", via Panisperna 89/A, 00184 Roma, Italy 



Contents 



1. Introduction 



1 



2. Higgs-pair production at the LHC 



2 



3. The decoupling limit of 2HDMs 



4 



4. Triple-Higgs self-couplings at the LHC 

4.1 Vector-vector fusion 

4.2 Higgs-strahlung 

4.3 Heavy-quark associate production 



7 
8 
10 



6 



5. Discussion of the results 



11 



6. Summary and outlook 



13 



1. Introduction 

The ability of the ATLAS and CMS experiments at the LHC to detect a SM Higgs boson 
over the full range of allowed masses has now been established by several detailed studies [1]. 
According to the latter, a multitude of production and decay channels will be accessible 
and will enable the determination of several of the couplings of the Higgs boson with 
accuracies that, depending on the precise value of the Higgs mass, Mh, can be as good 
as 10% [1]. Possible significant departures from the SM expectations will allow to infer 
that the underlying Higgs sector is not as simple as postulated in the SM. However, there 
is no guarantee that non-SM Higgs sectors will become manifest via these measurements. 
Examples are given by 2HDMs, including Supersymmetric models, where the spectrum of 
Higgs bosons beyond the lightest one could be very heavy and the couplings of the latter 
to fermions and gauge bosons reduce to those of the SM. We explore in this paper the 
possibility that in this limit a large deviation from the SM value X^ HH of the triple-Higgs 
self-coupling involving the lightest Higgs state, hereafter Xhhh, is allowed to survive. Such 
a circumstance would select the production of Higgs boson pairs as a possible channel for 
the identification of a non-SM Higgs structure. 

In Section 2 we review known results about Higgs-pair production in the SM. In Sec- 
tion 3 we discuss the main features of the decoupling limit in 2HDMs and demonstrate the 
existence of parameter sets leading to strongly enhanced triple-Higgs self-couplings, while 
maintaining the deviations of the couplings to fermions and to gauge bosons below the 
sensitivity range of direct LHC measurements. We verified that for (some of) these models 
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the analysis of Higgs-pair production can be performed rather reliably (for the channels 
we consider in the analysis) by simply rescaling in the Lagrangian the value of X HHH , 
keeping all other Higgs couplings fixed to their SM values 1 . Notice that this implies that 
our analysis is model independent 2 : our conclusions apply to any extended Higgs sector in 
regions of parameter space where the only observable departure from the SM is a sizeable 
deviation of Xhhh from X HHH • In Section 4 we focus on the study of Higgs-pair produc- 
tion at the LHC for IMH bosons, namely with 120 GeV < Mh < 140 GeV. In this region 
the detection of SM Higgs pairs is particularly challenging, due to the very large QCD 
backgrounds and the very small signal cross sections. We extend previous analyses by 
considering the production channels induced by vector-vector fusion, Higgs-strahlung and 
associated production with a top quark pair. For all these channels we assume the H^bb 
decay mode and we perform a background study. In Section 5 we present our assessment 
of the prospects to constrain deviations from the SM expectations, both at the LHC and 
at the so-called SLHC [2], the tenfold luminosity increase option of the LHC. Section 6 
contains our summary and conclusions. 

2. Higgs-pair production at the LHC 

The complete reconstruction of the Higgs potential necessarily requires the measurement of 
the Higgs self-couplings. These include a trilinear and a quartic interaction, parameterized 
by the coupling constants Xhhh and Xhhhh, which in the SM take the following values: 

A (o) __, M l (o) _oM H 

HHH ~ v ' HHHH ~ y 2 ' ^ > 

where v = 246 GeV is the Higgs doublet vacuum expectation value. A direct measurement 
of Xhhh can be obtained via the detection of Higgs boson pairs, wherein a contribution is 
expected from the production of a single off-shell Higgs boson splitting into two. However, 
the corresponding graphs will always be accompanied by diagrams where the two Higgses 
are radiated independently, with strength proportional to the Yukawa or gauge couplings. 
As a result, different production mechanisms will lead to different sensitivities of the Higgs- 
pair production rates on the value of Xhhh (notice that Xhhh contributes, through loop 
radiative corrections, also to vector-boson propagators and vertices [3], giving however 
effects that are too tiny to be detected even in precision observables). 

The leading production channels of Higgs boson pairs at hadron colliders [4] -[6] (see 
also [7]-[8]) are basically the same used for single-Higgs boson searches, namely: 

gg — »■ HH (gg fusion), gg,qq^QQHH (heavy-quark associated production), 



1 For such channels this is due to the fact that, if we isolate the diagram which involves the triple-Higgs 
self-coupling, its contribution is a constant at high energies. Therefore, despite the fact that Xhhh ^ Xhhh 
breaks gauge invariance, it is not crucial to enforce the unitarity of the cross section and it makes sense to 
assess the impact of changing the Xhhh value while keeping all other couplings fixed. We have explicitly 
verified, as discussed later, that this is indeed the case in the context of the 2HDM, namely, with a gauge 
invariant calculation. 

2 In this spirit, we will use the label H to identify a light Higgs state in both the SM and a generic 
extension, while reserving the symbol h specifically to the lightest Higgs boson of our example scenario, 
i.e., the 2HDM discussed in Section 3. 
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qqO — > qq^HH (vector-vector fusion), qq^—>VHH (Higgs-strahlung) , (2.2) 

wherein two Higgs bosons instead of one enter the final states and with V = or Z, 
Q = b,t and referring to any possible light (anti)quark flavor combinations 3 . The 
corresponding production rates are shown in Fig. 1 (borrowed from [7]), where we have su- 
perimposed those for heavy-quark associated production. The arrows indicate the variation 
in rate expected when changing the Higgs trilinear self-coupling in the range X^ HH /2 < 
Xhhh < 3/2 \$ HH - Depending on the value of Mh, different decay channels dominate the 
final state [8]. For a so-called IMH boson, with 120 GeV <M H < 140 GeV, the decay phe- 
nomenology is dominated by the channel H— > bb. For Mh > 140 GeV, the H—tW^^W^ 
and H—fZ^Z modes share the largest fractions of the decay rate. 

Given the rather low Higgs-pair production cross sections and the potentially large 
backgrounds associated with final states with the best decay rates (i.e., involving four 
6-quarks), naive arguments lead to the expectation that detection of IMH boson pairs 
via gg^HH^bbbb within the SM is most probably not feasible at the LHC and very 
difficult at the SLHC [2]. This had been suggested already in Refs. [9, 10] (see also 
Ref. [11]) and later quantitatively confirmed in Ref. [12], where the production and de- 
cay channel gg-^HH^bbT + T~ was shown to provide a better sensitivity to Xhhh- The 
potential of the rarer signatures gg^HH^bb^^) and gg^HH^bbfj, + fi~~ has been reviewed 
in Ref. [13], where the former was found to provide at the SLHC some limited scope in 
constraining Xhhh when Mh ~ 120 GeV. For a heavier Higgs boson, Mh > 150 GeV, 
the situation is much brighter in comparison. In Ref. [14] it was found that it should be 
possible at the LHC with design luminosity to establish that the SM Higgs boson has a 
non-zero trilinear self-coupling and that the ratio Xhhh/X hhh can ^ e restricted to the 
range 0-3.8 at 95% confidence level (CL), by exploiting gg^HH production and decay via 
HH^W ±( -*' ) W T W ± ^W T , in a variety of leptonic and/or hadronic final states. At the 
SLHC [2], such limits can be improved further and even 4-5<r excesses can be established. 

All such studies were based on the leading production channel of SM Higgs bosons, 
namely gluon-gluon fusion (see Fig. 1). We present here the first results of studies performed 
in the case of IMH boson pairs produced via the other three production modes in Eq. (2.2). 
The reason to exploit the latter is due to the additional triggers available in each case, with 
respect to the gg^HH mode: forward/backward jets in vector- vector fusion and high 
transverse- momentum leptons/light jets from gauge boson decays in both heavy-quark 
associated production and Higgs-strahlung. This could in principle enable one to perform 
additional independent measurements of Xhhh for an IMH boson, for which the current 
situation is extremely problematic. Although we will see that, as a result of our parton- 
level study, the available statistics is in general too low for quantitative estimates of X HHH , 
we will show that there is room for a signal, or for non-trivial limits, in the context of the 
2HDM discussed in the following Section. 

3 The case Q = b in heavy-quark associated production is actually irrelevant in the SM. In a general 
2HDM it may become significant. However, we will be interested only in the decoupling regime, which will 
be specified in more detail in the following Section. In this limit the departure of the cross section from the 
SM value is contrained to be at most within 70%. Hence, we have ignored the case Q = b altogether in our 
study. 
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Figure 1: Cross sections for Higgs-pair production in the SM via gluon-gluon fusion, vector- vector 
fusion and Higgs-strahlung (from [7]), plus top-quark associated production. The vertical arrows 
correspond to a variation of Xhhh from 1/2 to 3/2 of the SM value. 



3. The decoupling limit of 2HDMs 



A study of the decoupling limit of 2HDMs has recently been presented in [15] (see also 
[16]), where the general expressions for the spectrum and couplings of a generic, non-CP 
violating, 2HDM are derived. In particular, it is shown (see Eqs. (103)-(106) in [15]) that 
the approach to decoupling of the hVV, htt 4 and hhh vertices of the lightest 2HDM Higgs 
state h (of mass mh) can be parameterized as follows: 

(3.1) 
(3.2) 
(3.3) 



9hVV 
(0) 2 

9hvv 


~ 1 - 


AV 

4 ' 
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9htt 
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where the suffices H and (0) label the SM quantities, tua is the mass of the CP-odd Higgs 
boson, tan (3 is the ratio of expectation values of the up-type and down-type Higgs doublets, 



4 The expression for the case of 6-quarks is obtained from that for t-quarks upon the replacement cot f3 ■ 
tan (5. 
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A and A are dimensionless parameters of the 2HDM. Notice that the deviations from the 
decoupling limit are proportional to e = Xv 2 /m? A and e 2 in the case of the couplings 
to fermions and gauge bosons whereas the self-coupling is proportional to eA/A. The 
possibility that the ratio A/ A be large allows for the triple-Higgs self-coupling to remain 
large even when the other couplings are converging to their SM values. We analyzed this 
possibility by implementing the exact couplings of a generic 2HDM as discussed in [15] 
and scanning the parameter space in the range 1 < tan/3 < 50 5 , —Air < Aj < Air for all 
couplings Aj defined in [15], i = 1, ...7. While in the case A6,7 = full analytical constraints 
can be derived guaranteeing the Higgs potential be bounded [15], in the general scenario, 
A6,7 / 0, the stability of the potential has to be checked numerically for every point 
considered in the parameter space. Our general scan was subject to the constraints of 
tree- level unitarity [15] and to the requirement that the couplings g\yy, g\ tt and g\ hh differ 
from the SM values by no more than 30%, 30% and 70%, respectively. These values reflect 
the measurement accuracies expected after 300 fb _1 of accumulated LHC luminosity [l] 6 . 
The distribution of r = Xhhh/^hhh = ^hhh/^HHH f° r tne three Higgs mass values of 120, 
130 and 140 GeV in the general case (Ag,7 / 0) is shown in Fig. 2, where the scan assumed 
equiprobable input values for all 2HDM inputs. Some regions of the 2HDM parameter space 
can result in small values for run and itia- Hence, we have imposed a further constraint 
of no visibility at 3a level of the CP-even .ff-state at the LHC through the dominant gg 
fusion production channel, exploiting the projected significance of ATLAS and CMS for 
SM Higgs production. The gg^H production cross section has been estimated isolating 
the b- and t- loop contributions in the HIGLU [18] SM predictions and rescaling them with 
the proper 2HDM Higgs-quark couplings. As a result of this requirement, all the allowed 
points in Fig. 2 correspond to niH > 300 GeV and m,A > 250 GeV. (In the decoupling limit 
the dynamical behavior of the A state is very similar to that of H, so that one can safely 
assume that whenever the latter is undetectable the former is too.) From now on we shall 
refer to the region of parameters which survive the above constraints as decoupling limit of 
the 2HDM. The scan of all Aj leads to models with values of r in the ranges 7 

-8<r<36, M H = 120GeV, 

5 Note that in the scans we have traded the Higgs mixing angle a for the lightest Higgs boson mass, Mr, 
as independent 2HDM input parameter. 

6 Also a more optimistic scenario of 20%, 20% and 30% of measurement accuracies [17] has been inves- 
tigated, yielding conclusions similar to the ones outlined below, though over a restricted parameter range 
(see also next footnote). 

7 In the case of the more optimistic uncertainty scenario we get —3.5 < r < 18, —3 < r < 17 and —2 < r < 16 
for Mh = 120, 130 and 140 GeV, respectively. Notice that we did not include so far constraints on 
gl TT , because a full detector simulation for this case is still lacking, unlike the other cases. However, for 
completeness, we quote here some numbers also in presence of g\ Tr limits, though we will not adopt them 
in the remainder of the analysis. With a projected 10% accuracy on the latter, hence assuming 30, 30 
and 10% as constraints in the scan (with the error on g\ TT replacing the one on gi, ib ), the above numbers 
change to -1.5 < r < 18.6, -1.6 < r < 34.2 and -1.2 < r < 12.5 for M H = 120, 130 and 140 GeV, respectively. 
Assuming 20, 20 and 10% accuracies, the most optimistic case of all, we have —1.5 < r < 4.1, —1.6 < r < 12.4 
and — 1.2 <r < 12.1 for M H = 120, 130 and 140 GeV, respectively. Hence, as it will be clear from the 
following kinematical analysis, whichever the final outcome of the LHC analyses of single-Higgs production 
channels, all described scenarios for r will still be testable via double-Higgs production. 
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Figure 2: Distribution of the values of r = ^hhh/^hhh = ^hhh/^HHH m tne scans °f the 
2HDM parameters space defined in the text, for three values of the lightest Higgs boson mass. 
Normalization is to unity. 



-7<r<35, Mu = 130 GeV, 
-6 < r < 34, M H = 140 GeV. 



(3.4) 



4. Triple-Higgs self-couplings at the LHC 

This Section describes the general setup of our parton-level studies of signal and background 
rates. The numerical results are obtained by setting the renormalization and factorization 
scales to 2Mh for the vector-vector fusion and Higgs-strahlung signals and to V§ for heavy- 
quark associated production, with the latter choice accounting for the effect of the large top 
mass. For the various QCD backgrounds, the factorization/renormalization scale was the 
average jet Et for the production of light jets and \J!s for final states involving top quarks. 
Both Higgs processes and QCD backgrounds were estimated by using the parton distri- 
bution function (PDF) set MRST99(COR01) [19]. Both signal and background estimates 
were based on exact tree-level matrix element calculations using the ALPGEN program 
[20]. Furthermore, all signal rates have been cross-checked by using independent programs, 
based on either the HELAS subroutines [21] (for qq^^>qq^HH and qqC^^VHH) or the 
HELAC/PHEGAS package [22] (for gg^ttHH). As for numerical input values of SM pa- 
rameters, we adopt the ALPGEN defaults. We assume that 6-quark jets are distinguishable 
from light-quark and gluon ones. Unless otherwise stated, their tagging efficiency (includ- 
ing jet identification) is taken as = 50% for each b with E b T > 30 GeV and \rj b \ < 2.5 
(eb = otherwise) and we also require all light- and 6-jets to be separated, ARbb t bj,jj > 0.7. 
We do not assume fo-jet charge determination. Finally, in our parton-level analysis, we 
identify jets with partons. 
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Figure 3: Shapes of the integrated distributions of the minimum (left) and next-to- minimum (right) 
bb invariant masses in vector-vector fusion and relative QCD background, after the application of 
the cuts discussed in the text up to Eq. (4.1). Both signal and noise are shown for the case 
M H = 120 GeV. 



For the sake of definiteness, the signal results quoted in this Section will refer to the 
case of SM values for Xhhh- The analysis of the results in the context of our 2HDM 
example will be given in the next Section. 

4.1 Vector-vector fusion 

As seen in Fig. 1, within the SM, vector- vector fusion into Higgs boson pairs is the sub- 
process displaying the second largest cross section for any Higgs mass, behind gg^HH, 
and is the most sensitive to deviations of the Xhhh coupling from its SM value. In 
this channel, before any acceptance and selection cuts (and 6-tagging efficiency) one has 
production rates of about 1.4(0.7) [0.25] fb for M H = 120(130)[140] GeV. If one assumes 
efficient forward/backward jet tagging and high-purity sampling of 6-quarks, the dominant 
background is the QCD production of bbbbjj, where the two extra jets are required to have 
E 3 T > 20 GeV and to go one in the forward (rf > 2.5) and one in the backward < —2.5) 
direction. To enforce the reconstruction of the two Higgs bosons, we require all four 6's in 
the event to be tagged and that at least one out of the three possible double pairings of 
6-jets satisfies the following mass constraint: 

(m bub2 - M H f + {m b3M - Mh? < 2 a 2 m , (4.1) 

where a m = 0.12 Mjj- The bbbbjj cross sections with this event selection are 98(110) [130] fb 
for Mh = 120(130)[140] GeV (with no 6-tagging efficiency). The QCD background is dom- 
inated by configurations where the two 6-quark pairs are produced by the splitting of 
recoiling virtual gluons. These gluons tend to be off-shell by the smallest possible amount, 
compatibly with the requirement of AR separation and minimum Et of the 6-jets. The 
fake Higgs masses are most often reconstructed by pairing 6's travelling back-to-back in the 
transverse plane, namely pairs coming from the splitting of different gluons. As a result, we 
expect that even for events where two pairs exist passing the Higgs mass constraint, there 
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will be pairings of the four b's with rather low invariant mass. A very effective way to sup- 
press the QCD background is therefore to cut on the minimum and and next-to-minimum 
bb invariant masses, out of the six possible di-jet pairings, as can be appreciated from Fig. 3. 
The dependence of the signal and background cross sections on these additional constraints 
on rrifo masses, also including a more stringent cut on the Et of the forward/backward 
jets, are displayed in Tab. 1. For our forthcoming numerical investigation on the sensitiv- 
ity to anomalous triple-Higgs self-couplings we use the kinematic efficiency obtained with 
event selection c) of Tab. 1. The dependence of both signal and QCD background on the 
transverse energy of the forward/backward jets can be seen in Fig. 4, e.g., after the cuts 
in a) of Tab. 1. Other potential background processes, such as tibb, tiZ—>ttbb, ttH—^itbb, 
jjZZ^jjbbbb and jjbbZ^jjbbbb, have been found to altogether have cross section lower 
than 1CP 1 fb for event selection c) of Tab. 1 and have thus been neglected in the remainder 
of our study. The same is true for the process pp^bb + 4 light jets, where two light jets 
are mistagged as 6-jets with a fake efficiency of 0.01. 

Since in a realistic data analysis the background needs to be measured from the data, 
the question arises if the stringent cuts applied in the present analysis could produce 
peaking behaviours of the backgrounds in the signal region. To this aim we studied, for the 
most relevant background pp^Ab + 2 forward jets, the handle offered by the distribution 

-^r, with M = y/ (m ilti2 - M H y + (m i3 ,i4 - M H ) 2 (4.2) 

and {il,i2}, {i3,i4} the two pairings minimizing M. For the case Mjj = 120 GeV, the 
cut of Eq. (4.1) implies a cut M < 20 GeV. The distribution da/dM (see Fig. 5) does 
indeed exhibit a strongly peaked structure. However the peak of the distribution is outside 
the signal region at 120 GeV. We therefore expect that the background can be measured 
by looking at mass intervals away from the latter. As a cross-check of the accuracy of 
the extrapolation and as an assessment of the induced systematics, one can do the same 
exercise with a different nominal Higgs mass M, away from the signal region, for instance 
200 GeV (see also Fig. 5). This allows also the investigation of the region M < MV2a. 

Finally, concerning the trigger, in addition to the forward/backward jets, we have 
verified that more than 70% of the signal events are characterized by at least one fe-jet with 
Et > 100 GeV, which could give an additional handle in selecting Higgs-pair production 
via vector-vector fusion. 

4.2 Higgs-strahlung 

In the Higgs-strahlung channel, the cross section without any cuts and tagging efficiencies, 
before the decays of the and Z bosons, but including the HH^bbbb decay rate, is 
0.41(0. 19)[0.065] fb for M H = 120(130)[140] GeV. Let us first consider the case of leptonic 
decays of the gauge bosons, with t = e,/x. The cuts on 6-jets used here are the same 
as for the vector-vector fusion case, with the exclusion of the constraints discussed in 
Tab. 1. ZHH events are selected requiring two charged leptons with E T > 20 GeV, 
\rj \ < 2.5 and invariant mass within 5 GeV of the Z mass. W^HH events are extracted 
by isolating one charged lepton with E T > 20 GeV and \rf\ < 2.5 plus £vp lss > 20 GeV. The 
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Figure 4: Shapes of the differential (up) and integrated (down) distributions of the minimum 
forward/backward jet transverse momentum in vector- vector fusion and relative QCD background, 
after the application of the cuts discussed in the text up to Eq. (4.1) Both signal and noise are 
shown for the case Mh = 120 GeV. 
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Figure 5: Shape of the distribution of M in Eq. (4.2) for the case M H = 120 GeV (left) and 
M H = 200 GeV (right). 



resulting (cumulative) signal amounts to 0.0083(0.0048) [0.0019] fb, for M H = 120(130) [140] 
GeV, prior to 6-tagging efficiency. The corresponding irreducible background arises from 
Vbbbb events, with a rate of about 0.5 fb for the above cuts, with very small Higgs mass 
dependence. We also considered hadronic decays of both and Z, yielding a 46 plus 
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2(light-)jets signature. Acceptance, separation and mass cuts on the 46 system are as 
usual, whereas we require the two light jets to have E 3 T > 20 GeV, {rfl < 2.5 and to 
satisfy 70 GeV < rrijj < 100 GeV. The signal rates increase by a factor of about 4 with 
respect to the leptonic case, but the irreducible background increases by a factor of about 
10. Moreover, here, one should also include the large contribution coming from QCD 
production of jjbbbb, rendering the channel completely hopeless. 

4.3 Heavy-quark associate production 

As intimated already, in the SM only the 
case Q = t in Eq. (2.2) has some phe- 
nomenological importance. Taking our u- 
sual HH^bbbb signature leads to final sta- 
tes which include six 6-quarks, two of them 
coming from the decays of t and i. The 
largest backgrounds arise from the QCD 
production of ttbbbb and, if the request of 
tagging were limited to four 6-jets only, 
ttbbjj, where j represents a jet of gluons 
or light quarks. We estimated the ttbbjj 
background and obtained a rate of 1.5 pb, 
with usual jet acceptance cuts (Ej? > 20 
GeV, \r} b,3 \ < 2.5) and somewhat relaxed 
jet isolation requirements (AR^jjj > 0.4). 



M H (GeV) 


120 


130 


140 


Signal a) (fb) 


0.085 


0.049 


0.021 


Background a) (fb) 


40 


35 


32 


Signal b) (fb) 


0.073 


0.043 


0.018 


Background b) (fb) 


21 


18 


17 


Signal c) (fb) 


0.036 


0.021 


0.008 


Background c) (fb) 


4 


3.5 


3 



Table 1: The effects of additional cuts with re- 
spect to those described in the text up to Eq. (4.1) 
on the signal and background cross sections. Selec- 
tion a) means mj^ m > 50 GeV and m 



next— to — min 
bb 



> 



100 GeV, b) corresponds to mj£ m > 70 GeV and 
> 110 GeV and c) to the latter with 



rn 



next— to— min 
bb 



the more stringent cut E° T > 40 GeV (instead of 
E 3 T > 20 GeV) on the forward/backward jets. No 
6-tagging efficiency is included. 



This implies that tagging all six jets is mandatory. With the above cuts, the ttbbbb back- 
ground is in fact only 3.6 fb. Notice that, for this channel, maintaining a high 6-tagging 
efficiency is crucial, since it appears with the sixth power. If the SLHC luminosity led 
to a decrease from 50% to 40%, one would loose a factor of 4 in rate, and most of the 
advantage of the high luminosity gain would be lost. Notice also that in order to keep 
the ttbbjj background below the level of ttbbbb requires the fake tagging rate to be kept 
below 5% for each jet. To complete the reduction of the continuum background, we require 
that at least one out of all pairings of 6-quarks satisfies the constraint given in Eq. (4.1). 
We consider a semi-leptonic signature with one W ± decaying to leptons (e, n and their 
neutrinos) and the other one hadronically. So the final states are of the kind bbbbbbjj£i/£. 
Events are selected imposing E 3 T > 20 GeV, \rf\ < 2.5, ARj^jj > 0.4. We do not impose 
Et, i] or AR(j : ff, criteria on leptons from decay. The pre-selection criteria for 6-quarks 
are E\ > 30 GeV, \rf\ < 2.5, AR^ > 0.7. With the above cuts, we obtain the signal cross 
sections 0.021, 0.010 and 0.0035 fb (without including the branching ratio for the tt semi- 
leptonic signature, 8/27, but including the BR(H^bb) 2 ), corresponding to M# = 120, 
130 and 140 GeV, respectively. The irreducible background ttbbbb gives a cross section 
of 0.007 fb for any relevant Higgs mass. Both signal and background are expected to be 
sensitive to the choice of factorization/normalization scale, as they originate from QCD 
induced processes, primarily via gluon PDFs. Hence, we have checked the dependence of 
these numbers against variations of the factorization/normalization scale. For example, if 
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Q 2 = (2m t + 2M H ) 2 , we obtain instead the following signal cross sections: 0.025, 0.012 
and 0.004 fb for Mh = 120, 130 and 140 GeV, respectively, while the corresponding back- 
ground rate grows up to 0.013 fb. Detailed studies will be needed to assess to which extent 
a control data sample can be identified to fix the overall background normalization. The 
extraction of the signal from backgrounds could be feasible at a la statistical level, but the 
number of events is too small even to see a signal event, as the cross sections quoted above 
do not include the 66-tagging efficiency factor, i.e., ef. Assuming 56-tagging instead, always 
with €[> = 0.5, the signal events grow by a factor of 7, at the price of including further 
background processes: namely, tibb + light jets. With the usual fake tagging efficiency of 
1%, we checked that there is no sensible gain in statistical significance. 

5. Discussion of the results 

As was clear from the signal and background rates obtained in the previous Section, even 
the SLHC luminosity option will not be sufficient to establish a statistically significant 
signal for HH production within the SM in the 46 channel. In this Section we therefore 
consider the effects of an anomalous Xhhh coupling, in the ranges obtained in Section 3. 
In order to do so, we re-evaluated the signal rates in the presence of non-SM triple-Higgs 
self-couplings by simply rescaling the value of X^ HH to a generic Xhhh- The corresponding 
rates, as a function of Xhhh I^hhhi are snown m Fig. 6. Over the range allowed by our 
scan and by unitarity, cross section enhancements by up to two orders of magnitude can be 
obtained. Exclusion limits (at 95% CL) and signal evidence (at 3d) for anomalous triple- 
Higgs self-couplings by combining all channels are found in Tab. 2. In Tab. 3 the limits have 
been derived taking into account also a forward/backward jet recontruction efficiency of 
80%. The by far most sensitive mode is vector- vector fusion, with heavy-quark associated 
production being of some relevance only for small Higgs masses and positive Xhhh I X^ HH 
values, as can be expected by looking at Fig. 6. From the same plot, it is also obvious 
why the limits obtained via the Higgs-strahlung mode are always much weaker than those 
extracted in the other two channels. 

Notice that the described rescaling is not a gauge invariant operation, because diagrams 
involving other genuine 2HDM fields (namely, H, A and H^) in the Lagrangian, as well 
the rescaling of the Higgs-to-SM-particle couplings to the 2HDM values, are neglected. 
Therefore, one might suspect that the cross sections in Fig. 6 are anomalously enhanced 
by unitarity violation effects. To investigate this possibility we have re-computed the 
vector-vector fusion cross section, including the full set of 2HDM diagrams, for several 
hundreds of points uniformly distributed over the 2HDM parameter space, fulfilling the 
decoupling conditions of Section 3 and with r in the ranges accessible at the SLHC, namely 
r < -2.3(-3.1)[-5.3] or r > 3.9(5.9) [8.3], in correspondence of Mh = 120(130)[140] GeV. 

As a results of this exercise, we did find points, for any value of r, where the 2HDM 
cross section agrees almost exactly with the one depicted in Fig. 6. Since the calcula- 
tion is now gauge invariant, this fact alone demonstrates that the approximation implicit 
in Fig. 6 (namely neglecting additional graphs and coupling modifications) does not lead 
to an artificial enhancement of the sensitivity to Xhhh- Furthermore, we found a large 
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Figure 6: Dependence of the cross sections for the three processes qq( ^—>qq( ^HH, gg,qq-^ttHH 
and qcf ^—tVHH on r = Xhhh /^hhh> ^ n correspondence of three values of the Higgs mass, for 
the model setup discussed in the text. 



fraction (« 70%) of points where the cross section is reproduced by the above described 
approximation within a factor two. Finally, for a small, but sizeable, portion of points the 
cross section is substantially underestimated while for even fewer points is grossly overes- 
timated. For completeness, we also tried to understand the origin of such discrepancies. 
Firstly, the most noticeable ones are due to the BR(H^bb), as our definition of decoupling 
region allows for rather sizeable deviations of the latter from the SM value and, since we 
are looking for two Higgs bosons decaying into bb, this BR enters quadratically into our 
predictions. Secondly, also the diagrams proportional to gnw (V = W, Z) can give an 
important contribution: although the decoupling limit strongly constraints deviations of 
9hvv' from the SM value, the consequences of the latter can be enhanced by large destruc- 
tive interferences. Once these two effects are accounted for, the overall agreement is fairly 
good: the contribution of the SM diagrams to the total cross section is correct within a 
20% accuracy, except for a small fraction of points where the production rate is strongly 
underestimated. We have finally verified that the latter effect is due to neglecting diagrams 
where a light 2HDM Higgs pair is produced via the decay of the heavy 2HDM Higgs boson, 
H^hh. Notice that the effect of this additional contribution reinforces our conclusion that 
the signal is detectable, actually suggesting that a search is possible even in regions where 
Xhhh exhibits deviations from the SM value smaller than those appearing in Tabs. 2-3. 
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M H (GeV) 


120 


130 


140 


LHC, 95% CL 
SLHC, 95% CL 


-4.8 7.5 
-1.8 3.7 


-6.0 9.0 
-2.5 5.3 


-9.5 12.4 
-4.4 7.4 


LHC, 3a 
SLHC, 3a 


-6.6 9.3 
-2.7 5.1 


-8.1 11.0 
-3.6 6.5 


-12.4 15.4 
-5.9 8.9 



Table 2: Constraints on the ratio Xhhh /^ HHH using all the channels described in Section 4. 
(The results are almost completely driven by vector-vector fusion.) In the top box, the two values 
in each entry correspond to r m ; n , r max , where r < r m ; n and r > r max define the range which can 
be excluded at 95% CL (first row) or probed at the 3a level (second row), at both the LHC and 
SLHC). The number of signal events corresponding to 3a significance are about 120, 110 and 100 
for m,h = Mh = 120, 130 and 140 GeV, respectively, at the SLHC for the event selection c) of 
Tab. 1, with a background of about 1500, 1300 and 1100 events, respectively. 



M H (GeV) 


120 


130 


140 


LHC, 95% CL 
SLHC, 95% CL 


-5.6 8.0 
-2.2 3.8 


-6.9 9.8 
-3.0 5.7 


-10.8 13.7 
-5.1 8.0 


LHC, 3a 
SLHC, 3a 


-7.6 10.0 
-3.2 5.3 


-9.2 12.0 
-4.2 7.0 


-14.1 16.9 
-6.8 9.7 



Table 3: The same as in Tab. 2, including a 0.8 reconstruction efficiency for each for- 
ward/backward jet. 

In closing, we note that our results are not confined to the 2HDM in the decoupling 
limit but are truly model independent, thereby being applicable to other, more exotic Higgs 
sectors displaying a similar decoupling behavior between the lightest CP-even Higgs state 
and the heavier ones. The study of these models is however beyond the scope of this work. 

6. Summary and outlook 

In the SM, even at the SLHC, the measurement of X H hh for Mh < 140 GeV by using 
subleading Higgs-pair production channels, with the two Higgs particles both decaying 
into bb pairs, appears to be problematic, just like in the case of the leading mode studied 
in previous literature. However, by studying generic 2HDMs in the decoupling limit, we 
found large regions of parameter space where the triple-Higgs self-coupling Xhhh can 
differ considerably from the SM value, while keeping all other Higgs interactions with 
fermions and gauge bosons experimentally consistent with the SM limits. For the allowed 
regions, where the full 2HDM result can be emulated in a model-independent way (within 
an uncertainty of a factor of at worst two) by simply rescaling the SM triple-Higgs self- 
coupling to the 2HDM value, we have analysed the detectability (or otherwise) of various 
processes of the type pp^XHH^XAb in the mass range Mh < 140 GeV at the (S)LHC. In 
particular, signals at the level of 3a from non-trivial regions of the 2HDM parameter space 
can be obtained already at the LHC, with the SLHC extending the scope considerably 
further. The model-independent approximation used to perform the phenomenological 
analysis has been thoroughly tested against a complete 2HDM calculation for the case of 
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the vector-vector fusion process, which is the most sensitive to variations of the Xhhh 
coupling. The approximation can be improved (in a model-dependent way, however) by 
introducing the modified BR(H^bb) and coupling qhw (V = W, Z) as predicted in the 
2HDM. The only remaining discrepancies can be accounted for, within an accuracy of the 
order of 30%, by the contribution of the production of a heavy Higgs state and its decay into 
the light Higgs pair, being the contribution of the remaining non-SM diagrams negligible. 

In summary, Higgs-pair production at the (S)LHC could be an important channel to 
unravel a non-standard Higgs sector. In particular, in the mass window 120 < Mh < 140 GeV, 
the 4bX signature studied in this document will be the only accessible channel which could 
show sizeable departures of the trilinear coupling involving the lightest Higgs state from the 
SM value, in a scenario where the single-Higgs production channels do not show detectable 
anomalies with respect to the SM predictions. 

Finally, our further efforts will concentrate on the case of triple-Higgs self-couplings 
among neutral Higgs bosons in a general 2HDM away from the decoupling limit as well as 
in the Minimal Supersymmetric Standard Model (MSSM) (see [26] for 2HDM and MSSM 
studies of similar vertices involving charged Higgs bosons). 
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